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We provide a microscopic theory for the Doppler velocimetry of spin propagation in the presence 
of spatial inhomogeneity, driving electric field and the spin orbit coupling in semiconductor quan- 
tum wells in a wide range of temperature regime based on the kinetic spin Bloch equation. It is 
analytically shown that under an applied electric field, the spin density wave gains a time-dependent 
phase shift (j)(t). Without the spin-orbit coupling, the phase shift increases linearly with time and is 
equivalent to a normal Doppler shift in optical measurements. Due to the joint effect of spin-orbit 
coupling and the applied electric field, the phase shift behaviors differently at the early and the later 
stages. At the early stage, the phase shifts are the same with or without the spin-orbit coupling. 
While at the later stage, the phase shift deviates from the normal Doppler one when the spin-orbit 
coupling is present. The crossover time from the early normal Doppler behavior to the anoma- 
lous one at the later stage is inversely proportional to the spin diffusion coefficient, wave vector of 
the spin density wave and the spin-orbit coupling strength. In the high temperature regime, the 
crossover time becomes large as a result of the decreased spin diffusion coefficient. The analytical 
results capture all the quantitative features of the experimental results, while the full numerical 
calculations agree quantitatively well with the experimental data obtained from the Doppler ve- 
locimetry of spin propagation [Yang et ai, Nat. Phys. 8, 153 (2012)]. We further predict that the 
coherent spin precession, originally thought to be broken down at high temperature, is robust up 
to the room temperature for narrow quantum wells. We point out that one has to carry out the 
experiments longer to see the effect of the coherent spin precession at higher temperature due to 
the larger crossover time. 

PACS numbers: 72.25.-b, 72.25.Dc, 78.47.jj, 75.76.+j, 85.75.-d 



Understanding the spin transport phenomena is one of 
the most important issues in the fast developing field of 
spintronics [iHyl since it is crucial to the realization of the 
spintronic devices, such as spin-field-effect transistor 0- 
11 1 . In the proposed Datta-Das transistor [t|, the "on" 
and "off" states, distinguished by a 7r-phase difference in 
the spin precession mode, are switched by the gate volt- 
age which controls the coherent spin precession (CSP) of 
the passing carriers via the Rashba spin-orbit coupling 
(SOC) [12] acting as an effective magnetic field. To im- 
plement such devices, it is required that one is able to 
control and maintain spin polarization over a long enough 
distance, preferably at room temperature. 

Experimentally, the real space spin transport in semi- 
conductor is studied by using magneto-optic imaging 
13. 1 1-411 . or through conductance/current modulation 
lfj . 1 1 11 ] . An important development in the quantitative 
study of the spin transport is carried out by using spin 
transient grating spectroscopy [151 - 4201 ] . In these experi- 
ments, a spin density wave (SDW) with initial spin po- 
larization S z (x, 0) = So cos(qx) is created by two orthog- 
onal linear polarized light beams at time t — 0, where 
q is the wave vector of the SDW and x is the position. 
Under the influence of an applied in-plane electric field 
and the SOC, SDW picks up a phase and evolves into 
S z (x,0) = S°(q,t)cos(qx — <j)(t)). By optically monitor- 
ing the temporal evolution of the amplitude S® (q, t), one 



obtains the spin diffusion coefficient and relaxation rate 
l3-21|. Very recently, the spin drifting and CSP were 
studied by the Doppler velocimetry which monitors the 
phase shift 4>{t) [15|,l22j. For pure spin drift without spin 
precession, the phase shift is simply <p(t) = qv c it, where 
Vd stands for the drift velocity. The linear increase of 
phase shift with time is equivalent to a Doppler shift 
Aoj = Vdq. When both the CSP and drifting are present, 
the phase shift deviates from this simple relation and be- 
haves anomalously. In Ref. [TH Yang et al. reported that 
at low temperature (T — 30 K), <j>(t) deviates strongly 
from the simple relation of qvat and clearly shows the 
anomalous behavior caused by the CSP. However, once 
the temperature rises to 150 K, the anomalous behavior 
of the phase shift disappears in the time frame of the ob- 
servation. It was then concluded that CSP breaks down 
at high temperature, although the mechanism of the dis- 
appearance of CSP is not clear 15 1. 

The temporal evolution of the SDW contains all of the 
important information of spin transport, such as spin dif- 
fusion coefficient, spin mobility and CSP. The transient 
spin grating spectroscopy together with the Doppler ve- 
locimetry therefore enable one to quantitatively study the 
spin transport in semiconductors. However, the dynam- 
ics of SDW is quite complex when spatial inhomogeneity, 
applied electric field as well as the SOC are present. To 
correctly extract the information from the experimen- 
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tal data, a thorough understanding of spin transport is 
in demand. Without the electric field, the amplitude 
S®(q,t) decays double exponentially when the SDW dif- 
fuses along the [110] crystal axis in (001) GaAs quantum 
wells as spin rotates along the net effective magnetic field 
due to the SOC and diffusion 0, [H |23|. Based on the 
kinetic spin Bloch equation (KSBE) approach [f| [23 - |29| , 
it is shown that the information about spin diffusion co- 
efficient, CSP and spin relaxation can be extracted from 
the wave- vector dependence of the two decay rates 21 1. 
In this paper we will further extend the theory to include 
the electric field and show, both analytically and numeri- 
cally, that the theoretical and experimental results agree 
well with each other. We will demonstrate that the CSP 
is robust even in high temperature regime for narrow 
quantum wells and point out that to observe the effect 
of the CSP at high temperature, one has to carry out 
the observation for a longer time than the case at low 
temperature. 

As will be shown later, the SDW transporting along 
the [110] crystal axis with the wave vector q in a (001) 
GaAs quantum well evolves as 



S z (x,t) 



S z {q,0) expKA? 2 + l/r s )i]/2 
X [e-^^ cosiqx-Vdiq + q^t] 

+ e 2D ™ t cos[qx-v d (q-q Q )t]} , (1) 

when the SOC is weak enough. Here D, t s and Vd are 
the spin diffusion coefficient, spin relaxation time of the 
SDW and the drift velocity under the electric field E, 
respectively, qo = m*(j3 + a) and q' = m*(p + a) with 
m* representing the effective mass, a being the Rashba 
coefficient [l2| and j3, $' both standing for the coeffi- 
cients of the linear Dresselhaus term [30j with correc- 
tions from the cubic Dresselhaus terms [21(. Without 
the applied electric field, Vd — and the amplitude of 
the SDW decays double exponentially with fast and slow 
rates Dq 2 + l/r s ± 2Dqq [l3,l2l|. When there is an ap- 
plied electric field but without the CSP (qo = q'o = 0), 
the SDW decays exponentially and gains a phase shift 
which changes linearly with time with a slope qvd- In 
experiments, the phase shift can be monitored optically 
and the slope qvd is equivalent to a normal Doppler shift. 
With the CSP, the situation is more complex. For small 
time t <C 1 / (4:D\qqo\) , the fast and the slow modes share 
the same weights. Therefore the phase shift of the SDW 
reads 

(j)(t) ~ [v d (q + q' )t + V d (q - q' Q )t}/2 = v d qt . (2) 

Nevertheless, for large time t ^> l/(4_D|<7<7o|), the slow 
mode dominates and the phase shift becomes 



Doppler one without the CSP. After some time, the slope 
deviates from the normal one and reduces to v d (q — q'o)- 
The sign of the slope reverses when q < q' . In the special 
case of q = q' , the phase approaches a stationary value 
at large time. The crossover time of (j)(t) from the normal 
Doppler behavior in the early stage to the anomalous one 
at the later stage is about t c ~ \/{AD\qqo\). In the case 
of small diffusion coefficient, wave vector of SDW or the 
SOC, this crossover time can be larger. 

Equation (p} captures all of the qualitative features of 
the experimental results [lj| at low temperature, from 
the two modes in the temporal evolution of SDW when 
E = (hence v d = 0), to the details of how the phase 
changes with time and wave vector when E ^ 0. Specifi- 
cally, the temporal evolution of the phase can be divided 
into two stages. In the early stage, the phase increases 
with time with a steeper slope. The slope becomes flatter 
in the later stage, or even reverses its sign for small wave 
vectors. The larger the wave vector is, the quicker the 
phase behavior changes from the early stage to the later 
stage. All these features qualitatively agree with Eq. ([T]) . 

As a result of the increa sing electr on-p honon scatter 



ing and spin Coulomb drag [lj|2j,|3lj-|33|, the spin diffu- 



~ v d (q - %)t 



(3) 



That is, in the presence of the CSP, the phase shift first 
changes linearly with a slope qvd, same as the normal 



sion coefficient D decreases with the increase of temper- 
ature. It is therefore expected that at high temperature 
the crossover time is larger than that at low temperature. 
With the correction of the cubic Dresselhaus term, q also 
reduces as the temperature rises and the crossover time 
is further prolonged. If the contribution from the cubic 
Dresselhaus term is so large that qo approaches zero, then 
the CSP breaks down completely at high temperature. In 
the recent experiments, the CSP is thought to be broken 
down at high temperature based on the lack of the clear 
anomalous behavior in the phase shift from the observa- 
tion in a limited time regime. The cubic Dresselhaus term 
is speculated to be the main cause for the breaking down 
[lij . However, this is highly unlikely because qo at low 
temperature only differs by a few percents from the one 
at high temperature for the narrow quantum well used in 
the experiments. It should be pointed out that, since the 
behaviors of phase in the early stage are the same with or 
without the CSP, one should be cautious in determining 
the existence of such CSP from the experimental data in 
limited time regime, especially when the crossover time 
t c is large. In our opinion, purely from the existing exper- 
iments, it is inconclusive to determine whether the CSP 
survives or not at high temperature. Moreover the ana- 
lytical solution presented in this paper and the previous 
numerical calculations support that CSP is stable even 
at room temperature [2l|, [2J-[29| . 

To derive the solution [Eq. ([T])] and to determine if the 
CSP is stable at high temperature, we turn to the full 
KSBEs for the spin transport in a (001) GaAs quantum 
well of well width a grown along the z-axis @, @, HI , 24 - 
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FIG. 1: (Color online) Phase shifts of drifting spin grat- 
ing under an applied electric field of E = 2 V/cm for dif- 
ferent wave vector at (a) T = 30 K and (b) T = 150 K. 
Blue Curve/Triangle: q — 1.07 cm -1 ; Green Curve/Square: 
q — 0.81 cm" 1 ; Red Curve/Diamond: q — 0.59 cm -1 ; Purple 
Curve/Circle: q = 0.34 cm -1 , respectively. Inset of (a): Fit- 
ting of spin relaxation times r+ (Red Curve/Circle) and r_ 
(Blue Curve/Triangle) at T = 30 K. The dashed purple curve 
in (b) is the phase shift for q = 0.34 cm" 1 at T = 300 K. 
All the curves in the figures are from theoretical calculation 
whereas the symbols are the experimental data from Ref. HH . 



8p(x, k, t) 



eE(x) 



dpk(x,t) k x dp k (x,t) 



dt Bk x 
i[h tot (k))-<r/2,^ k (a:,i)] 



to* dx 
dp k (x,t) 



Of 



(4) 



Here we assume that the transport direction is along the 
x-axis. pk(x,t) are the density matrices whose diago- 
nal elements fu.<r(x,t) represent the electron distribution 
functions with spin a and momentum k = (k x ,k y ) = 



(kcoscj),ksm(j)) at position x. The off-diagonal elements 
of Pk(x, t) stand for the spin correlations between spin-up 
and -down electrons. Equation includes drift of elec- 
trons driven by the electric field E(x) (the second term), 
diffusion caused by the spatial inhomogeneity (the third 
term), spin precession around the total magnetic field 
ht t(k) (the fourth term), as well as all the scattering 
(the last term). The total magnetic field is composed 
of the external magnetic field B in the Vogit configura- 
tion, the effective magnetic field h(k) due to the Dressel- 
haus/Rashba SOC as well as the one from the Hartree- 
Fock term of the electron-electron Coulomb interaction. 
The expressions for the Hartree-Fock and the scattering 
terms are given in detail in Refs. [jj l34l [35| . 

By expanding the density matrix p{x, k, t) = 
J2 l pi(x,k,t)e u '^, the KSBEs can be expanded as a se- 
ries of coupled equations for pi(x,k,t). By neglecting 
the Hatree-Fock term 37| and the inelastic scattering 
and using the fact that the spin density is S(x,t) = 
Tr{crp(a;, k, t)} = J °° Tr[crpo(x, k, t)]kdk/2TT, one 
finds that, to the leading order, under a uniform applied 
electric field E the spin density obeys the following equa- 
tion 



OS _ D (PS 
dt dx 2 



- BS BS 
2Dh x — + v d — + v d ti xS-f' S,(5) 
ox ox 



in which the diffusion coefficient D = (k 2 ''t\ / '(2m* 2 )) , the 
drift velocity Vd = (eEri/m*). For transport along the 
[110] crystal axis (cc-axis), as in the experimental setup, 
h = —%Gy and h' = —q' e y with e y being the unit vector 
along the y-axis ([110] crystal axis), qo = m*(/3 + a) 
and q' = m*0' + a). Here (} = - ~j{k 2 )/A and 
P' = f3 — ry{k 2 )/2 with fj — 771- 2 /a 2 and 7 being the Dres- 
selhaus coefficient [3(J. Noted that we have included the 
corrections from the cubic Dresselhaus term. The relax- 
ation matrix S% = d±ag{l/T Xl 1/t v , 1/t z }, with 1/t x — 



((/3+a) 2 fc 2 r 1 /2) + (7 2 fc 6 T 3 /32), l/r y = (^- a fk 2 r 1 /2) + 
(7 2 fc 6 T3/32) and 1/t z = 1/t x + 1/t v being spin relaxation 
rates of spin components along the x, y and z direction, 
respectively. 1/tj = J Q 27r T ^ g - ) cos(W)d9 /2-k with r(k,9) 
standing for the momentum relaxation time due to the 
electron-impurity scattering. For a system near the equi- 
librium, (•••) = / •••a/(e k )/9e k d 2 k/ / Bf( £]i )/Be k d 2 k 
with /(e) being the Fermi distribution function. 

The first two terms of the right hand side of Eq. ([5]) 
describe the spin diffusion and CSP around the net ef- 
fective magnetic field, proportional to h due to the joint 
effect of the SOC and diffusion respectively. The third 
and fourth terms show the two effects of the electric field 
on the dynamics of the SDW, i.e., drift of electrons under 
the influence of driving electric field and the CSP around 
another net effective magnetic field, Vdh. 1 , due to the joint 
effect of the SOC and drifting. The last term describes 
the spin relaxation of the SDW. It is noted that a simi- 
lar equation is derived by the random walk approach at 
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zero temperature [20|. However, the diffusion coefficient 
differs from the one obtained in this paper by a factor of 
2, and the correction from the cubic Dresselhaus term is 
only partially taken into account in Ref. [20j . 

For SDW with wave vector q, one can simply replace 
d x by iq in Eq. ([5]) and obtains the equation for S(q,t), 
the Fourier component of SDW. For initially z polarized 
SDW, the solution reads 

S z ( q ,t) = S z ( q ,Q)^ + (Q)e- r+{q)t + X-(q)e- r - {,})t }, (6) 

with S z (g, 0) standing for the initial spin density, 

r±(«) =Dq 2 -iv d q + l/T s ±A/2T y , (7) 

\±{q) = [l±l/A]/2, where 1/r. - (1/r, + l/r 2 )/2 and 
A = yjl + ATjj(2Dqqo — ivdq'^) 2 - The result is similar 

to the one without the applied electric field: the tem- 
poral evolution of S z (q,t) is composed of two modes, a 
fast decay mode with decay rate K{r + (q)} and a slow 
one with rate 5ft{r_(g)} [2l|. For small electric field, 
the change in the decay rates is quadratic to the field. 
Therefore, the electric field has only marginal effect on 
the decay rates. However, the electric field introduces ad- 
ditional spin precession, causing S z (q,t) to oscillate with 
frequency |3{r±(g)}| which is linear to the electric field 
for small field. In the case with |A| ;§> 1, such as in the 
system with weak SOC or in the special case when /3 ~ a, 
the solution can be further simplified to A± ~ 1/2 and 

r±(g) = Dq 2 ± 2Dqq - iv d (q ± q' ) + l/r s , (8) 

one then gets time evolution of the SDW as Eq. (p}. 

To clarify whether the CSP survives at high tempera- 
ture where qo reduces and electrons suffer strong inelastic 
scattering, in addition to the elastic scattering, one has 
to calculate the temporal evolution of the SDW by nu- 
merically solving the full KSBEs. In the calculation, we 
include all the relevant scattering, such as the electron- 
impurity, electron-AC phonon, electron-LO pho non as 
well as electron-electron Coulomb scattering [38]. The 
material and structural parameters are chosen from the 
available experimental data directly or by fitting the cor- 
responding experimental data: quantum well width a = 
9 nm, and the electron density N e = 1.9 x 10 11 cm -2 . The 
effective impurity density iVj is set to be 0.12N e by fitting 
the mobility under the laser intensity of 0.25 /zJ-cm~ 2 
[39| | . The Rashba coefficient a = since the quantum 
well is symmetric. The Dresselhaus coefficient 7 is set 
to be 12 meV-A -3 by fitting the spin relaxation time at 
T = 30 K. The KSBEs used here are valid unless the 
higher subbands of quantum well are significantly occu- 
pied by electrons. For the parameters used here, this 
does not happen until T ~ 800 K, well above the room 
temperature. For wave vector-dependent spin relaxation 
times at T = 30 K, shown in the inset of Fig. HJa), the 
theoretical and experimental results are also in very good 



agreement. To see the effect of the CSP, we plot the 
phase shifts <p(t) of the SDW under the influence of an 
applied electric field of E — 2 V/cm at different gratings 
as function of time in Fig. Q]for T = 30 and 150 K. For 
comparison, we also plot the experimental data in the 
figure. It can be seen that theoretical and experimen- 
tal results agree well with each other. More importantly, 
at T = 150 K both theoretical and experimental results 
show influence of the CSP on <fi near t — 150 ps, where 
the slopes of <f> become flat. In the theoretical calcu- 
lation, the effect is revealed more clearly for SDW with 
q = 0.34 x 10 4 cm -2 after about 200 ps when the crossover 
from positive slope to the negative one is completed. The 
calculations at room temperature is also carried out and 
CSP is found to be robust even at room temperature. 
But the crossover from the early to later stage is fur- 
ther delayed, e.g., for q = 0.34 x 10 4 cm -1 the crossover 
from positive slope to the negative one finishs at around 
250 ps, as shown in Fig. [ljb). From these calculations 
one concludes that CSP indeed survives at room tem- 
perature for the narrow quantum wells studied here, but 
one has to carry out the observation for a longer time to 
observe its effect on the phase experimentally. 

In conclusion, we provide a microscopic theory for the 
Doppler velocimetry of spin propagation in the presence 
of spatial inhomogeneity, driving electric field and the 
CSP due to the joint effect of the SOC and transport in 
a wide range of temperature regime. Applying this the- 
ory to study the transport of SDW, we analytically show 
that in the presence of the electric field the SDW gains a 
time-dependent phase shift <t>(t). Without the CSP, <fr(t) 
grows linearly with time with slope qvd for the SDW with 
wave vector q, which is equivalent to a normal Doppler 
shift in optical measurements. Due to the CSP caused 
by the net effective magnetic field from the joint effect of 
the SOC and transport, the short time and the long time 
behaviors of phase shift are different. At the early stage 
<j){t) grows with time in the same way as that of without 
the CSP, i.e. with a normal Doppler slope qvd- At the 
later stage the slope reduces to (q — qo)vd, deviating from 
the normal Doppler one. For small q the slope at the later 
stage even reverses its sign. The crossover time from the 
early to the later stages is inversely proportional to the 
spin diffusion coefficient D, the wave vector of the SDW 
q and the SOC strength. Since D decreases with the 
increase of temperature, the crossover from the early to 
the later stage at high temperature would be prolonged. 
Our numerical calculations, which include all the rele- 
vant scattering such as the electron-impurity, electron- 
AC phonon, electron-LO phonon and electron-electron 
Coulomb scattering, agree quantitatively well with the 
existing experimental results. By extending the calcu- 
lation time beyond the experiment regime, we predicate 
that the CSP, originally thought to be broken at high 
temperature, is robust and stable up to the room tem- 
perature. We further point out that to observe the effect 
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of the CSP on the phase shift a longer observation time 
is required at higher temperature. 
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